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HWhy  Quantum  Chemical  Reaction 

Kinetics  Studies? _ 

■  Only  Option  When  Experiments  are  not  Possible  or  Limited 

■  Combustion  Conditions  of  P  &  T  too  Extreme  to  Probe 

■  Accuracy  (Ea)  can  be  as  Good  or  Better  Than  Experiments 

■  Thermochemical  Accuracy  Possible 

■  Ideal  for  Branching  Ratio  Predictions  for  Closely  Competing  Reactions 

■  Can  be  a  Cost  Effective  Alternate  to  Experiments 

■  Hardware  &  Software  Efficiencies  Improving  Constantly 

■  A  Balanced  Approach  to  Kinetics  Calculations  Recommended 


Microscopic  Rate  Coefficient 


Electron  Correlation 


Zador  et  al.,  . 
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Our  Interest  in  Combustion 

Chemistry 


■  Autoignition  Chemistry 

■  Low  Temperature/Low  Pressure  Conditions 
■  Bi-propellants 


■  Decomposition  Chemistry 

■  Emerging  Energetic  Materials 

■  Explosives 

■  Ionic  Liquid  Propellants 
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N  H  +  NO  Reaction  Kinetics 


Perspective 


Radical  Chemistry  Modelling 

■  N2H4/N02  Autoignition 


■  Recent  Works 

■  Only  Theoretical  Studies 

■  See  Raghunath  etal.,  Adv.  Quantum  Chem.,  69,  253  (2014) . ^298K,  iatm=  (2.3  x  10'11) 

■  See  Daimon  et  a/.,  Sci.  Tech.  Energetic  Materials,  74,  1 7  (2013) . k298  K  1atm  =  (1 .6  x  1 0'14)  ? 

■  See  Daimon  etal.,  J.  Propul.  Power,  30,  707  (2014) . k298  k, iatm=  (1-9x  10-11) 

■  Also,  See  Kanno  etal.,  DOI:  10.1021/acs.jpca.5b00987  (2015) 

CH3NNH2  +  N02  -»•  Products,  k  =  (2.2  x  10'13) 
trans-CH3NHNH  +  N02  -►  Products  k  =  (1.4  x  10'12) 
cis-CH3NHNH  +  N02  — >  Products  k  =  (1.2  x  10'12)  cm3  molecule'1  s_1 

■  This  work 

■  Pulsed  Laser  Photolysis  -  Flow  Tube  MS  Experiments 

■  First  Experimental  Determination 

■  Ab  initio  Chemical  Kinetics 

■  Multi-reference  Second-order  Perturbation  and  Coupled-cluster  Methods:  PES  (Potential 
Energy  Surface) 

■  RRKM  (Rice-Ramsperger-Kassel-Marcus)  Theory  and  Master  Equation  Simulations:  k 
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Pulsed  Laser  Photolysis 
Flow  Tube  Apparatus 


Pump 


QMS, El 


DAS 


Flow-tube  reactor 


Vacuum  chamber 


Gas  mixture 
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N2H3  Source  & 
Flow  Tube  Chemistry 


N2H4  +  hv  ^N2H3  +  H 
N02  +  H  NO  +  OH 
N2H4  +  H  N2H3  +  H2 
N2H4  +  OH  -►  N2H3  +  H20 
N2H3  +  N02  -»•  N2H2  +  HONO 
N2H3  +  N02  — ►  other  products 
HONO  -»■  loss 
N2H3  — >  loss 
H  — ►  loss 
OH  — ►  loss 


0193  nm  =  450  x  IO-20  cm2  molec1 
k2  =  1.3  x  IO10  cm3  molecule1  s_1 
k3  =  1.4  x  10'13  cm3  molecule1  s'1 
k4  =  3.6  x  10'11  cm3  molecule1  s_1 

k5 

k6 

k7=  1  s-1 
k8  =  10  s-1 
k9  =  10  s-1 
k10=  10  s-1 
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[HONO]  =  (k5[N02][N2H3]0){e  k7  ‘  -  e  *’ *}/(*’  -  k7) 
k’  =  (k5  +  k6)[N02]  +  k8 


[N2H4]  =  5  x  io14,  [H]  =  5  x  1012,  [N02]  =  1  x  io13  to  5  x  1013  molecule  cm  3 
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signal  (arb  units) 


n2h3  +  no2 

Reaction  Kinetics 


Typical  [HONO]  Temporal  Profile  Second-order  Plot 


k298 K, 2 Torr n9  =  (1-23  ±  0.25)  x  IQ-11  cm3  molecule'1  s-1 
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Potential  Energy  Surface  for 
NoHoNOo  Adduct  Formation 


CASPT2/CBS 
RCCSD(T)/CBS//CASPT2 
RCCSD(T)/CBS//B3LYP 
Unit:  kcal/mol 
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Potential  Energy  Surface  for 
NoHoONO  Adduct  Formation 


a 
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Branching  Rate  Coefficients  (2  Abmr 
No)  &  High  Pressure  Limit  Versus  T 
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NH2NHN02,  high  pressure  limit 
NH2NHN02 

trans-N2H2+trans-HONO 

cis-N2H2+transHONO 
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MNB  &  DNB  Energetic  Materials 

■  MNB & DNB 

■  Mononitrobiuret  &  1 ,5-dinitrobiuret  (less  stable) 

■  promising  explosives 

■  Recent  Works 


■  See  Geith  et  al.,  Propellants,  Explosives,  Pyrotechnics,  29,  3  (2004) 

■  AHcomb(DNB)  =  (51 95  ±  300)  kJ  kg-1  (bomb  calorimetry  and  MP2/cc-pVTZ  A Hf)  cf  HMX  9435  &  RDX  9560  kJ  kg-1 

■  Vd  =  8660  ms-1,  cf  HMX  9100  &  RDX  8750  ms'1 

■  See  Geith  et  al.,  Combust  and  Flame,  139,  358  (2004) 

■  Recent  synthesis  (known  since  1898  by  Thiele)  &  decomposition  mechanism  studies 

■  Initially  NH2N02  released,  followed  by  residues  decomposing  to  give  HNCO 

■  See  Liu  etal.,  J.  Phys.  Chem.  A115,  8064  (201 1 )  &  Sun  etal.,  J.  Phys.  Chem.  A118,  2228 
(2013) 

■  Quasi-classical  direct  dynamics  trajectory  simulations  to  understand  DNB  decomposition 

■  Elimination  of  HNN(0)OH  intermediate  identified 

■  See  Suntsova  et  al.,  Struct.  Chem.,  24,  745  (2013) 

■  Electronic  structure  of  DNB  was  studied  by  two-dimensional  B3LYP  potential  energy  scans 

■  Solid-  and  gas-phase  conformers  differ 


■  This  work 

■  PESs  for  the  thermal  decomposition  of  MNB  and  DNB  investigated 

■  Temperature  and  pressure-dependent  rate  coefficients  were  calculated  using  micro-canonical 
transition  state  theory  with  master  equation  simulations 
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MNB  &  DNB  Structures 
(M06-2X/aua-cc-pVTZ  level! 


MNB  (trans)  MNB  (cis) 


H-bond  Stabilized  Conformers 

■  MNB(CS)  0 

■  MNB(trans)  +4.82 

■  MNB(C1)  +6.38 

■  MNB(cis)  +7.95  kcal/mol 


■  Double  H-bond  Stabilized  Conformers 

■  DNB(CS)  0 

■  DNB(C2v)  +1.22 

■  DNB(C-i)  +7.82  kcal/mol 


DNB  (Cs) 


DNB  (C-,) 


DNB  (C2v) 
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Stationary  Point  Energies  for  MNB 
Thermal  Decomposition 


HBC-P1 


IM1  +  HISINIOIOH 


HBC-P2 


I  M2  +  HNCO 


trans-biuret  radical  +  NO, 


HBC-P3 


IM3  +  HNCO 


NH-.C(0]NH,+  MO.NCO 


iso-trans-biuret  radical  +  NO, 


0.00 

33.86 

18.50 

29.42 
34.58 
27.75 

30.43 

51.41 

6.41 
41.05 
34.23 
40.33 
53.77 
34.38 
57.25 
7.84 


0.00 

33.18 

20.51 

26.85 

33.73 

26.20 

32.17 

51.34 
6.38 

40.35 
32.65 
38.46 
54.71 
32.53 
56.99 
8.29 


0.00 

33.99 

19.43 
27.08 
34.66 
26.41 
32.95 
52.75 
5.86 
41.00 
32.47 
38.77 
53.71 
30.68 

58.43 
7.32 


All  energies  include  zero-point  corrections  and  are  in  kcal/mol  relative  to  the  global  energy 
minimum  of  MNB  (Cs). 


Within  1.0  kcal/mol  for  the  barrier  heights,  and  1.8  kcal/mol  for  endothermicities  compared  to 
RCCSD(T)/cc-pV®Z//M06-2X/aug-cc-pVTZ,  implying  that  the  M06-2X/aug-cc-pVTZ  energies 
are  good  for  larger  analogous  systems. 
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CT||^  Potential  Energy  Surface  for  MNB 
■UsM  Decomposition 

RCCSD(T)/cc-pVooZ//M06-2X/aug-cc-pVTZ 
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Ab  initio  Kinetics  of  MNB  Thermal 

Decomposition 


Terrperature  1DDG/T(K) 

High-pressure  rate  limit  for  MNB  (Cs)  and  MNB  (Ct) 
primary  thermai  decomposition  channels 


A,  -  =  t .32  x  io13T° 131  exp  (-17474.7 fT) 
k2„=  1.05  x  io13T0112  exp  (-17692.1/T) 
kg  od  =  1.15  x  1014T^17exp  (-18011.4/T) 
k4  .  =  1 .26  x  104  T2  539  exp  (-21 375.5/T) 
As'  -  =  2.07  x  1023T-2  002  exp  (-27714.8/T) 


MNB  {CJ  and  MNB  (C1>  primary  thermal 
decomposition  channels 


Below  1200  Kt  HNCO  &  HNN(Q}OH  eliminations 
dominate,  consistent  with  experimental  work  of 
Gieth  eta/,2004 

Above  1200  Kr  N02  elimination  dominates 


10* 


0  5  10  15  2.0 

Temperatire  I000/T(K) 


Rate  coefficients  at  pressure  of  1  atm  for  MNB  {CJ 
and  MNB  (C,)  primary  thermai  decomposition 
channels 


Pressure  dependence  of  unimofecular  decomposition 
evaluated  using  AE^^  =  200  *  (T/300)0  &5  cm-1  energy 
transfer  probability  model 

Strong  pressure  dependences  observed 

HNCO  &  HNN{0)OH  eliminations  are  competitive  & 
dominate  below  2500  K 

N02  elimination  becomes  competitive  only  above  2500  K 

R2  &  R3  {HNCO]  &  Rt  {HNN(O)OH)  are  primary  (product) 
channels  at  ignition  temperatures 
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Stationary  Point  Energies  for  DNB 
Thermal  Decomposition 


Species 

M06-2X/aug-cc- 

pVDZ 

M  06-  2X /a  ug-cc-p  VIZ 

DNB  |CJ 

0.00 

0.00 

TS6 

34.17 

33.77 

HBC-P6 

21.58 

20.30 

HNN{0)0H  +  N02NHC,(0)NCO 

28.68 

26.79 

TS7 

31.72 

31.21 

HBC-P7 

30.70 

29.98 

IM2  +  N02NCG 

34.79 

36.91 

Cj-MNB  radical  +  N02 

50.20 

49.79 

trans-MNB  radical  +  N02 

50.89 

51.09 

DNB  (CJ 

1.16 

1.22 

TS8 

28.64 

28.61 

HBC-P8 

22.35 

21.70 

HNN{0)0H  +  iso-NOjN  HC(0  )NCO 

29.08 

27.26 

cis-MNB  radical  +  N02 

52.54 

52.84 

DNB  (C,) 

7.45 

7.82 

TS9 

9.42 

10/32 

TS10 

10.79 

11.20 

All  energies  include  zero-point  corrections  and  are  in  kcai/mot  relative  to  the  global  energy  minimum  of  DNB  (Cs). 

Predicted  energy  barrier  heights  from  these  two  levels  of  calculation  show  excellent  agreement,  but  discrepancy 
reaches  to  2.1  kcal/mol  for  reaction  endothermicities.  Consequently,  the  M06-2X/aug-cc-pVTZ  energies  were 
used  for  carrying  out  the  kinetics  analysis  for  DNB  decomposition. 
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Potential  Energy  Surface  for  DNB 

Decomposition 


M06-2X/aug-cc-pVTZ 
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Ab  initio 


Kinetics  of  DNB  Thermal 
Decomposition 
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Tempe  nature  101 

High-pressure  rate  limit  for  DNB  ( 
primary  thermal  decompositi 


k6  _  =  3.22  x  10wT-° 239  exp 
k7  „  =  8.49  x  10l6T-° 632  exf 
kg  „  =  9.72  x  1 013  T-°  248  exp 
kg  .  =  9.96  x  1023  T-2205  expWi  ^ou.o,  ■ , 


0.5  1.0  1.5  2.0 

Temperature  1000/T(K) 

Rate  coefficients  at  pressure  of  1  atm  for  DNB  (C 
primary  thermal  decomposition  channels 


DNB  (TS7) 


ig  pressure  dependences  observed 
w  1200  Kr  HNN(0)OH  elimination  dominates 


-4 


VefJ  - 577 


MNB(TSI) 


DNB  (TS8) 


1,574 


MNB  (TS2) 


2vt 


Below  1500  Kt  HNN(0)QH  elimination  dominates 
from  DNB  (C^ 

Above  1500  K,  NG2  elimination  dominates 
from  DNB  (C2v) 

Similar  trend  seen  in  Liu  etai.,  2011,  and  good 
agreement  above  1500  K 


k3  for  HNN{0}OH  elimination  from  DNB  is  much  larger 
than  kt  for  HNN{0)OH  from  MNB  below  800  K 

Consistent  with  TS1  (and  TS3)  energy  higher  than  TS8 
{Also,  TS2  higher  than  TS7) 

TS3  is  stabilized  by  additional  intramolecular  H-bondlng 
{Also,  TS7  is  a  late  transition  state) 


DNB  is  less  stable  than  MNB 

Present  quantitative  interpretations 
are  consistent  with  experimental  work 
of  Gieth  ef  at.  2004 
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Conclusions  & 
Acknowledgements 


■  N2H3  +  N02  addition  reaction  is  fast 

■  trans-N2H2  +  trans-HONO  main  products 

■  PES  surface  characterized 

■  Calculated  reaction  rate  coefficients  in  agreement  with  experiments 

■  MNB  &  DNB  stabilized  via  6-member-ring  moieties  involving  intramolecular  H-bonding 

■  Solid  state  conformers  have  smallest  dipole  moment 

■  Energy  barriers  and  endothermicities  at  the  M06-2X/aug-cc-pVTZ  level  of  theory  show 
remarkable  agreement  with  the  values  obtained  from  RCCSD(T)/cc-pV°°Z//M06-2X/aug-cc- 
pVTZ  computations 

■  Former  level  of  theory  should  be  applicable  to  larger  analogous  systems 

■  MNB  decomposition  initiated  by  the  elimination  of  HNCO  and  HNN(0)0H,  the  latter  is  also 
released  in  DNB  (C2v)  decomposition 

■  Energy  barrier  for  HNN(0)0H  elimination  in  DNB  is  6.60  kcal/mol  lower  than  that  in  MNB 
due  to  an  extra  hydrogen  bond  in  the  transition  state  for  the  former 


■  DNB  less  stable  than  MNB  (as  previously  shown  by  experiments) 


